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ABSTRACT 

The effect of velocity gradient on stagnation region heat transfer 
augmentation by free stream turbulence was investigated. Heat 
transfer was measured in the stagnation region of four models 
with elliptical leading edges with ratios of major to minor axes of 
1:1, 1.5:1, 2.25:1, and 3:1. Four geometrically similar, square bar, 
square mesh, biplane grids were used to generate free stream 
turbulence with different intensities and length scales. Heat 
transfer measurements were made for the following ranges of 
parameters: Reynolds number, based on leading edge diameter, 
37,000 to 228,000; dimensionless leading edge velocity gradient, 
1.20 to 1.80; turbulence intensity, 1.1 to 15.9%; and length scale 
to leading edge diameter ratio, 0.05 to 0.30. Stagnation point heat 
transfer augmentation by free stream turbulence can be predicted 
using a modified version of a previously developed correlation for 
a circular leading edge. Heat transfer augmentation was 
independent of body shape at the stagnation point The heat 
transfer distribution downstream from the stagnation point can be 
predicted using the normalized laminar heat transfer distribution. 

INTRODUCTION 

The heat transfer distribution around a streamlined object 
immersed in a flow usually has a maximum in the stagnation 
region. Cooling the stagnation region is important in many 
industrial application; however, none is more critical than in the 
modem high efficiency gas turbine. Combustor exit temperatures 
often exceed the melting temperature of superalloy turbine vane 
materials. Accurate prediction of stagnation region heat transfer 
is vital to turbine designers. 

For a laminar free stream flow, heat transfer in the stagnation 
region can be predicted using Frossling’s solution [1] if the 


pressure distribution around the object is known from say an 
invidd calculation. Unfortunately, in a turbine as in most other 
internal flows, the flow is not laminar. Combustor primary and 
dilution jets and the wakes from upstream blades and vanes cause 
high levels of turbulence. A turbulence intensity of 11% was 
measured at the exit of a combustor [2]; modem high enthalpy 
rise combustors probably produce even higher levels. Free stream 
turbulence can augment stagnation region heat transfer, ratios of 
turbulent to laminar stagnation heat transfer of 1.9 have been 
measured [3]. 

Stagnation region heat transfer augmentation in the presence of 
free stream turbulence is believed to be caused by vortiaty 
amplification (see [4] for a review). Free stream turbulence can 
be viewed as a continuum of tangled, vortical filaments. 
Filaments that are convected into the stagnation region, with 
components normal to the stagnation line and normal to the free 
stream flow, are stretched and tilted by the divergence of 
streamlines and acceleration around the bluff body. This 
stretching causes the vortiaty to be intensified through 
conservation of angular momentum. It has been shown both 
experimentally and numerically [5,6,7] that voitidty in the 
stagnation region causes heat transfer to be increased while the 
boundary layer remains laminar. 

It has been known for many years that free stream turbulence 
can augment stagnation region heat transfer [8,9]; however, 
results of experiments are inconsistent and attempts to correlate 
heat transfer augmentation as a function of turbulence intensity 
and Reynolds number, while ignoring the length scale, 
[10,11,12,13,14] have not been entirely successful. 
Any resulting correlations usually predict the author’s data but not 
data from other researchers. 

Lowery and Vachon [15] measured lateral length scale in 
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their study of the effect of grid generated turbulence on stagnation 
region heat transfer but they did not have a sufficient variety of 
grids to deduce an effect of scale. Their resulting correlation has 
been used as a standard against which subsequent data sets have 
been compared, sometimes with large discrepancies; see for 
example [16]. 

There have been several attempts to isolate the effect of 
turbulence length scale; Yardi and Sukhatme [17] used four 
different grids to generate a range of length scales. The four grids 
were all of different geometry, i.e. two were screens and two were 
triplane grids, all had different rod spacing to rod diameter ratios. 
They showed a trend of increasing heat transfer with decreasing 
length scale; however, there was so much scatter in the data that 
their daim of ten boundary layer thicknesses for an optimum 
length scale is questionable. 

Dyban et al [18] used perforated plates as well as a fully 
developed turbulent pipe flow to investigate the effect of intensity 
and scale on stagnation region heat transfer. Their results showed 
increasing augmentation with decreasing scale but they did not 
attempt to correlate the data based on this finding. 

More recently, Ames [19] used simulated combustor segments 
to generate turbulence and measure its effect on heat transfer to 
a flat plate and a stagnation region. Ames concentrated on 
relatively large scale turbulence where the length scale to leading 
edge diameter ratio was greater than 1.0. He used the rapid 
distortion theory of Hunt [20] and the measurements of Hunt 
and Graham [21] near a plane surface to develop a new 
correlating parameter involving Reynolds number, turbulence 
intensity, and what Ames calls an energy scale (the average size 
of the energy containing eddies). Ames used three different 
diameter cylinders to investigate stagnation region heat transfer, 
his data were correlated well using his new parameter. The data 
of several other researchers were also correlated by his parameter 
but with more scatter. 

Van Fossen and Ching [22] measured heat transfer on a 
circular leading edge downstream of four different square bar, 
square mesh turbulence generating grids. They developed a 
correlation for stagnation region heat transfer involving turbulence 
intensity, integral length scale, and Reynolds number that fit their 
results to ±4% and the results of other authors with similar 
turbulence generators to within ±8%. 

Considering vortitity amplification theory, it would seem 
reasonable that leading edge velocity gradient would have an 
effect on stagnation heat transfer. Higher velocity gradients 
would cause more rapid stretching of the vortical filaments as 
they are convected past the leading edge thus increasing heat 
transfer. The purpose of the present work was to study the effect 
of leading edge velocity gradient on stagnation region heat 
transfer augmentation. Three models with elliptical leading edges 
were fabricated with heat transfer gages in the stagnation region. 
The ratio of major to minor axes for the models were: 1.5:1, 
2.25:1 and 3:1. All three models had the same leading edge 
radius of curvature as the circular leading edge of Van Fossen and 
Ching [22]. Each of the models was qualified in a low turbulence 
flow by comparing stagnation heat transfer measurements with a 
numerical solution for laminar stagnation flow. The same four 
grids used in [22] were used to generate turbulence upstream of 


each model. The grids were square mesh, biplane grids made 
from square bars with different bar widths. Each of the four had 
identical mesh to bar width ratio. 

Stagnation region heat transfer was measured with each grid at 
various distances upstream of each of the models. Data were 
taken at Reynolds numbers based on leading edge diameter 
ranging from 37,000 to 228,000. Turbulence intensities were in 
the range 1.1 to 15.9 percent while the ratio of integral length 
scale to leading edge diameter ranged from 0.05 to 030. 
Stagnation point velocity gradient varied from 1.20 to 1.80. 

It will be shown that stagnation heat transfer augmentation due 
to turbulence is unaffected by the velocity gradient near the 
leading edge and can be predicted for each of the elliptical bodies 
by the same correlation developed for the circular leading edge 
[22]. A method for determining the heat transfer distribution 
downstream of the stagnation point will also be presented. 

NOMENCLATURE 

A surface area of gage, m 2 

a constant in equation (9) 

a e ellipse major axis 

B ratio of model thickness to tunnel height 

b turbulence grid bar width, cm 

b e ellipse minor axis, cm 

C constant equation (11) 

constant in equation (7) 

C T parameter in equation (1) 

d diameter of model leading edge = 2 R, cm 

E mean hot wire voltage, V 

Cfws fluctuating output voltage from linearizer, V 

Fr(s/R) Frossling number 
/ constant in equation (10) 

k air thermal conductivity, W/mK 

M mesh spacing of bars in turbulence grid, cm 

m exponent in equation (9) 

n exponent in equation (7) 

Nu Nusselt number 

Pr Prandtl number 

q heat flow, W 

R leading edge radius, cm 

Re Reynolds number 

Rf x) autocorrelation of velocity signal 

r recovery factor 

s surface distance from stagnation, cm 

T temperature, °C 

Tu turbulence intensity 

U mean velocity, m/s 

u’ st ream wise RMS fluctuating velocity component 

x st ream wise distance, cm 

Greek symbols 

A integral length scale, cm 

p air density, Kgm/m 3 

a standard deviation 

x time shift, s 


2 


<D 


heat transfer augmentation factor, equation (12) 
Subscripts 

0 evaluated at zero flow 

avg average 

b bar width 

d leading edge diameter 

El electrical heating 

gap epoxy filled gap between gages 

lam laminar free stream 

r recovery 

rad radiation 

st static 

/ total 

turb turbulent free stream 

w wall 

x streamwise 

oo free stream 


TEST FACILITY & INSTRUMENTATION 
Wind Tunnel 

The experiments were carried out in the wind tunnel shown in 
Fig. 1, which is described in detail by Van Fossen and Simoneau 
[6]. Air drawn from the test cell passes through a flow 
conditioning section and a 4.85:1 contraction before entering the 
15.2 cm wide by 68.6 cm high test section. The maximum 
velocity attainable was about 46 m/s ec. Gear tunnel turbulence 
levels were less than 0.5 percent for all flow rates. After leaving 
the test section, the flow passed through a transition section into 
a 10-inch pipe which contained a standard ASME orifice run with 
flange taps. Air then passed through a butterfly valve which was 
used to control the tunnel flow rate and then to the Laboratory 
exhaust system. The readings from four Chromel-Constantan 
thermocouples located around the perimeter of the inlet to the 
conditioning section were averaged to yield the stagnation 
temperature. 

Turbulence Grids 

For the present tests, four square bar, square mesh, biplane 
turbulence generating grids were used. The grids were fabricated 
keeping the ratio of mesh spacing to bar width constant at 45 
yielding an open area of 605 percent Grid parameters are 
defined in Fig. 2 and dimensions of the grids are given in Table 
I. Henceforth grids will be referred to by the symbol given in the 
table. Turbulence generating grids could be installed at axial 
locations ranging from 2.41 to 52.3 cm upstream of the model 
stagnation point allowing length scale and intensity to be varied. 

Heat Transfer Models 

The four heat transfer models used in this study had elliptical 
leading edges. The ratio of major to minor axes, ajb c , were 1:1, 
1.5:1, 2.25:1 and 3:1. All models had the same radius of 
curvature, R, of 330 cm at the stagnation point A comparison 
of the model profiles is shown in Fig. 3. All models had wedge 


shaped afterbodies that extended about 61 cm downstream of the 
leading edge to eliminate vortex shedding. Fig. 4 is a photograph 
of the heat transfer models and afterbodies. The purpose of the 
4 models was to provide different velocity gradients in the 
stagnation region to determine if this would have an effect on 
stagnation heat transfer augmentation. Leading edge velocity 
gradients calculated using an invicid 2-D panel code [23] are 
shown in Fig. 5. Velocity gradients made dimensionless by 
leading edge radius and free stream velocity ranged from 1.2 to 
1 . 8 . 

The circular model had nineteen heat flux gages and all the 
elliptical models had twenty-nine heat flux gages embedded 
symmetrically around the stagnation line. Fig. 6 is a sketch of a 
typical model cross section showing the heat flux gage 
arrangement Each beat flux gage consisted of an aluminum strip 
6.60 cm long by 0.476 cm wide and 032 cm deep. A Kapton® 
encapsulated, foil, electric heater was fastened to the back of each 
aluminum strip with pressure-sensitive adhesive. The temperature 
of each gage was measured by a Chromel-Alumel thermocouple 
embedded in a groove. A guard heater behind the heat flux gage 
array prevented heat conduction to the interior of the model. The 
average gap between the aluminum strips was 0.025 cm and was 
filled with epoxy. The dimensionless surface distance, s/R, from 
the stagnation line to the center of each gage is given in Table II. 
The aluminum strips were maintained at a uniform constant 
temperature by a specially designed control circuit, (see Van 
Fossen et al., [24]). Steady state, spanwise-averaged heat 
transfer coefficients were calculated for each aluminum strip 
based on the power supplied to the strip and the wall-to-fluid 
temperature difference. 

Instrumentation and Data Acquisition 

Turbulence intensities, and integral length scale were measured 
using a standard, 5 pm, single, hot-wire oriented perpendicular to 
the flow direction. A spectrum analyzer computed 
autocorrelations which were then used to determine the length 
scale. 

Steady-state operating conditions (temperatures, pressures, 
voltage and current to gages, eta) were recorded on the 
Laboratory data acquisition system called ESCORT [25]. For 
every heat transfer data point, twenty readings of each data 
channel were recorded. These twenty readings were averaged to 
give a single value for each channel. To eliminate any offset 
between data channels, a reading was obtained by shorting all the 
inputs to ESCORT and subtracting this "zero" from each 
subsequent reading. 

EXPERIMENTAL PROCEDURE 
Turbulence Parameters 

Hot wires used in the turbulence measurements were calibrated 
in an open air jet at the same temperature as the wind tunnel 
flow. Velocity calibrations used a two point, iteration method in 
conjunction with the analog linearizers as described in [26]. 
The system frequency response was estimated to be around 30 
kHz with the standard square wave test. Turbulence intensity and 
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integral length scale were measured at several locations 
downstream of each of the four grids and at several tunnel 
velocities without the heat transfer models present. 

Heat Flux Measurements 

For the heat transfer measurements, all the heat flux gages were 
heated to temperatures of about 46 °Q the average recovery 
temperature of the air was around 27°C giving wall to air 
temperature difference of approximately 19°C. This temperature 
difference was a compromise; large enough to keep errors in 
temperature difference small and small enough to minimize 
thermal property variations. All of the heat flux gages were 
maintained at the same temperature within ±0.1°C; this was 
accomplished by adjusting the gain of each control circuit to 
maintain the gage thermocouple voltage within one microvolt 
Heat flux measurements were carried out with each grid installed 
at several axial locations from the stagnation line of the leading 
edge. For each grid position, tests were performed at three 
Reynolds numbers ranging from 37,000 to 228,000. 

DATA REDUCTION 
Turbulence 

Turbulence intensity downstream of the grids was calculated as 
the ratio of the linearized root mean square (RMS) of the 
fluctuating component of hot wire bridge voltage to the mean 
voltage. 

Integral length scale was determined from the autocorrelation 
function, R(t), by fitting it with an exponential function using 
least squares 

R(x) = e <x W 


Fr{s!R) = (<?g — 

A(T w - T r )k 


where q B is the heat added by the heater (voltage x current), q ^ 
is the heat lost by radiation, and q^ is the heat conducted away 
to the epoxy gap and to the unguarded ends of the gages. A is the 
exposed heat transfer gage surface area, T w is the gage 
temperature, T r is the recovery temperature at the gage location, 
and lb is the thermal conductivity of air. 

An estimate of the gap loss, q^, can be obtained from an exact 
solution for two-dimensional heat conduction in a rectangle half 
the epoxy gap width wide and the aluminum gage depth deep. 
Two adjacent sides are assumed insulated, one side held at the 
constant temperature of the aluminum strip and the final side 
convecting to the air at the local recovery temperature. Heat 
conducted out the unguarded ends of the gages can be estimated 
from the same analysis by assuming a large gap width. Details 
of this analysis are given in [24]. 

Corrections for radiation heat loss, were made assuming 
gray body radiation to black surroundings and an emissivity of 
0.05 for the aluminum gage. Heat lost through the sides and 
unguarded ends of the strips was on the order of 10 percent of the 
total heat flow, while the radiation heat losses were on the order 
of 0.2 percent. 

The recovery temperature was calculated from 

t, = r„. ♦ <r, - (4) 


Data between 0.33 s; R(x) £ 1.0 were used for the curve fit 
Integrating equation (1) over time delays ranging from 0 to oo, 
and using Taylor’s hypothesis that time delay and streamwise 
distance are related by the mean velocity, U y the integral length 
scale then becomes 


a , - u£s(T)d x = JL ( 2 ) 


where T^ m is the static temperature upstream of the model and 
T t is the total temperature. The recovery factor, r, was calculated 
as [27] 


r = 1 - 


p m 

(p£0. 

v v '/ 


(1-y/pT) 


( 5 ) 


This method eliminated the problem of determining the upper 
limit of integration for autocorrelation functions that oscillate 
about zero for large time delays. More details of the turbulence 
measurements are presented in [27]. 

Heat Transfer 

Power from the electric heaters is removed from the aluminum 
strips by convection to the air, radiation to the surroundings, and 
conduction to the epoxy gap between the gages where it is 
convected to the air. The Frossling number was determined from 
an energy balance on each gage 


where the mass flow ratio, pU(s)/(pU) m was found from a 
numerical solution of the flow over the model; the solution 
included the effects of the tunnel walls [28]. 

The thermal conductivity, viscosity, and Prandtl number were 
evaluated at the free stream total temperature from equations 
given in [29]. Total temperature was used to evaluate the 
thermal properties because in [28] a numerical study showed that, 
if the thermal properties were based on a reference temperature 
that involved the wall temperature, reversing the direction of heat 
flux (cooling the wall) caused an undesirable change in the 
Frossling number. 

The Reynolds number. Re# was based on the diameter of the 
leading edge, d , and the mass-velocity averaged between the flow 
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area with maximum model blockage and the unblocked upstream 
flow area, i.e. 


(p£0 * 


= (p v). 


(2 ~B) 
2(1 -B) 


( 6 ) 


where B is the ratio of maximum model thickness to tunnel height 
and ranged from 0.096 for the circular leading edge to 0.293 for 
the 3:1 ellipse. 


UNCERTAINTY ANALYSIS 


RESULTS AND DISCUSSION 


Turbulence 

Turbulence intensity was measured as a function of distance 
downstream of each grid without the heat transfer model in place. 
The data was in general agreement with that of Baines and 
Peterson [31]. Each grid and Reynolds number had slightly 
different characteristics so intensity data for each case were fit 
with a power law curve of the form 



(?) 


Twenty samples were obtained for each steady-state 
measurement and averaged to minimize random errors. Standard 
deviations were also obtained from the twenty samples and used 
as an estimate of random error. Estimates of the accuracy of each 
measuring instrument were then made, added to the random 
component, and combined by the method of Kline & McQintock 
[30]. Results of the uncertainty analysis indicated an average 
uncertainty of ±65% for the Frossling number. The contributions 
from individual measurements to the overall uncertainty in the 
Frossling number are shown in Fig. 7. 

Error in turbulence intensity was estimated by assuming that the 
linearizer approximates King’s law, i.e. the velocity could be 
expressed in terms of bridge voltage as 



k c *~ j 


where E is the bridge output voltage, E 0 is the voltage at zero 
velocity, and n and are constants. Differentiating this 
expression and dividing by the velocity, one obtains an expression 
for turbulence intensity 


u' m dU ^ 2nEe r^ 
~U m ~ E 2 - El 


( 8 ) 


where dU is taken as the RMS of the fluctuating component of 
velocity, and dE has been replaced by the RMS of the 
fluctuating component of bridge voltage, e^s- The exponent, n, 
was assumed to be near 2.0 with an error of ±10%. The method 
of error estimation described above was then applied to this 
expression; typical uncertainties estimated for the turbulence 
intensity measurements were on the order of ±15%. Uncertainties 
of the turbulence length scale were assumed to be the same order 
as the turbulence intensity. 


Coefficients for each of the curve fits appear in Table HI. 
Turbulence intensity was found to vary by less than 5% in the 
spanwise direction. 

In [22], X-wire measurements are reported that showed the 
turbulence from the square bar, square mesh grids to be nearly 
isotropic for x/b greater than about 25. With a few exceptions for 
grid Gl, all the beat transfer data was obtained with the models 
greater than 25 bar widths downstream of the grids. 

Roach [32] developed a correlation for the integral length 
scale of grid generated turbulence of the form 


A fx\ (10) 

b - nIW 


Length scale data from the present grids were found to have the 
same square root of distance dependence but the coefficient, /, 
varied from grid to grid and was an average of 35% larger than 
the value found by Roach. The coefficients for the length scale 
correlation are also found in Table III. 

When the model is present downstream of the grid, turbulence 
is distorted as the stagnation point is approached. The fluctuating 
component of velocity increases and the mean velocity approaches 
zero [6] sending intensity levels very high. This brings up the 
problem of where to evaluate the turbulence intensity and length 
scale for use in a heat transfer correlation. It was felt that tests 
of most turbulence producing components, e.g. combustor, would 
be conducted without the model present; therefore, turbulence 
intensity and length scale used in the following correlations were 
evaluated from the curve fits in Table III using the distance from 
the grid to the stagnation point of the model. Turbulence 
intensity varied from 1.1 to 15.9%. The ratio of length scale to 
leading edge diameter ranged from 0.05 to 0.30. 

Heat Transfer 

Verification. Heat transfer results in the leading edge region 
with no turbulence grid in the tunnel are shown in Fig. 8 for the 
four different models. Measured frees tream turbulence intensity 
in this case was less than 0.5 percent. The ordinate for the heat 
transfer plots is the Frossling number, Nu^Re^) 172 . Data are 
presented as a function of surface distance from stagnation made 
dimensionless by the leading edge radius, R. In all cases the data 
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agree, to within the estimated experimental error, with the 
Frossling solution obtained using velocities calculated from a 
panel code [23] and with a 2-dimensional numerical solution from 
the PARC-2D code [28] thus confirming the accuracy of the 
experimental technique. The worst agreement between the 
experimental and numerical results is for the 2.25:1 model, where 
the experimental results are from 1.4 to 9.2% above the numerical 
results at the stagnation point. The upper limit of discrepancy is 
above the estimated experimental error. This model seemed to 
have a mind of its own; some days the model gave results that 
agreed quite well with the numerical results and other days large 
errors were observed. Possible causes of this error were 
investigated including model profile, surface irregularities, 
thermocouple calibration, and dear tunnel turbulence level. The 
model was X-rayed to see if the internal guard heater was 
touching the surface heat flux gages. No obvious cause could be 
found; therefore, the data for this model is presented "as is". 

Stagnation Point Augmentation. In [22], the stagnation heat 
transfer data for the drcular leading edge was correlated by the 
function 


Fr{ 0) = 0.00799 



/ \ 

-0574 

rr n 0*800 

Tu Re d 

$ 

1 





+ C 


(ii) 


The constant, C, was set to the zero turbulence Frossling number 
of 0.939 which was determined from the PARC-2D calculation. 
The other constants were determined from a least square fit of the 
data. The function was found to correlate the data to within ±4%. 

In [28], it was shown by numerical calculation that the 
stagnation point Frossling number with simulated turbulence 
(sinusoidal veloaty variation upstream of the leading edge) 
divided by the laminar Frossling number was independent of body 
shape. Following this line, the correlation for the circular leading 
edge was modified by dividing by the laminar stagnation Frossling 
number. The modified correlation then gives the stagnation point 
heat transfer augmentation factor, <£, due to free stream turbulence 




= 0.00851 

~ 0 0500 

Tu Re d 

/ \ 
\ 

-0574 

% 





( 12 ) 
+ 1 


Values for the terms FiiO)^ for the four models are 0.939, 0.870, 

0.811, 0.775 in order from the circular leading edge (1:1 ellipse) 
to the 3:1 ellipse; these values were also taken from the PARC-2D 
numerical solutions. 

Comparison of the correlation for stagnation heat transfer 
augmentation by free stream turbulence (equation (12)) and the 
experimental data is shown in Fig. 9. The correlation was 
developed using only the circular leading edge data, thus, the fit 
for that data is the best. In general, the fit is excellent falling 


mostly within the ±4% bands drawn on the figure. The 2.25:1 
model has the most scatter, as mentioned earlier, this model had 
problems. If the Frossling number for the 2.25:1 model had been 
normalized using the average of the experimental low turbulence 
data instead of the numerical solution, agreement with equation 
(12) would have been much better. The excellent agreement for 
the other three models confirms the validity of this correlation 
method. 

Equation (12) contains no term that involves the velocity 
gradient at the stagnation point; yet the stagnation point heat 
transfer augmentation is the same for several different levels of 
velocity gradient Thus, the hypothesis that heat transfer 
augmentation above laminar levels should increase in the presence 
of higher velocity gradients is not bom out by these experiments. 


Distribution of Heat Transfer Around Leading Edge. Fig. 
10 is a plot of the local Frossling number normalized by the 
stagnation value versus dimensionless surface distance from the 
stagnation point for each of the models. The symbols represent 
an average of the local Frossling number data for all turbulent 
free stream conditions (all grids, Reynolds numbers, and grid 
positions). The dotted lines represent the standard deviation of 
the normalized data and the solid line is the PARC-2D solution 
for a laminar free stream which has been similarly normalized. 
Agreement between the normalized turbulent heat transfer 
distribution and the normalized laminar distribution is good; thus, 
a good prediction of the heat transfer at a given distance from the 
stagnation point can be obtained by using the correlation 
developed to predict the stagnation heat transfer and multiplying 
by the ratio of local to stagnation heat transfer from a solution for 
the laminar free stream. 


CONCLUSIONS 

Spanwise average stagnation region heat transfer measurements 
have been made on four models with elliptical leading edges 
downstream of turbulence generating grids. The ratio of major to 
minor axes for the elliptical leading edges ranged from 1:1 to 3:1; 
all the models had the same leading edge radius of curvature. 
Velocity gradients at the stagnation point ranged from 1.20 to 
1.80. Four turbulence generators were used, they were square 
mesh, square bar, biplane grids with identical mesh spacing to bar 
width ratios and bar widths ranging from 0.16 to 1.27 cm. 
Reynolds numbers based on leading edge diameter ranged from 
37,000 to 228,000, turbulence intensities ranged from 1.1 to 
15.9%, and the ratio of integral length scale to leading edge 
diameter ranged from 0.05 to 030. Results are summarized as 
follows: 

1. The Frossling and PARC-2D solution and the experimental 
stagnation region heat transfer data for a laminar free stream are 
in good agreement for all models tested. 

2. Stagnation point heat transfer augmentation by free stream 
turbulence for the elliptical leading edges can be predicted using 
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a normalized version of a previously developed correlation for 
a circular leading edge. 

3. Dimensionless heat transfer augmentation due to turbulence 
is independent of body shape and therefore independent of 
velocity gradient at the stagnation point 

4. The heat transfer distribution downstream from the 
stagnation point can be predicted using the normalized laminar 
heat transfer distribution and the stagnation point heat transfer 
correlation. 
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Table I. Turbulence generating grid dimensions. 


GRID NO. 

b, cm (in) 

M/b 

% OPEN AREA 

G1 

1.270( .500) 

4.5 

60.5 

G2 

0 • 635 ( .250 ) 

4.5 

60.5 

G3 

0.318( .125) 

4.5 

60.5 

G4 

0 . 159 ( .063) 

4.5 

60.5 


Table II. Heat flux gage dimensionless surface distances from stagnation point. 


MODEL 

U1 1 

LliiiU 

| 2 .25 1 1 | 

1 3>1 

GAGE 

s/R 

1 

0.000 

0.000 

0.000 

0.000 

2 

0.152 

0.151 

0.152 

0.151 

3 

0.303 

0.303 

0.304 

0.301 

4 

0.455 

0.456 

0.457 

0.452 

5 

0.607 

0.607 

0.610 

0.604 

6 

0.759 

0.759 

0.762 

0.755 

7 

0.910 

0.910 

0.915 

0.907 

8 

1.062 

1.062 

1.068 

1.058 

9 

1.214 

1.214 

1.221 

1.210 

10 

1.365 

1.366 

1.374 

1.361 

11 


1.518 

1.527 

1.512 

12 


1.670 

1.680 

1.664 

13 


1.822 

1.833 

1.816 

14 


1.975 

1.986 

1.968 

15 


2.127 

2.140 

2.119 
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Table III. 


Power law curve fits of turbulence intensity and integral length scale data. 


- - * nr 


K 

b 



Grid 

Velocity 

symbol 

Re b 

a 

m 

I 

G1 

R1 

38650 

206.1 

-0.875 

0.240 

G1 

R2 

18000 

206.1 

-0.875 

0.240 

G1 

R3 

7934 

206.1 

-0.875 

0.240 

G2 

R1 

17190 

146.3 

-0.780 

0.272 

G2 

R2 

9514 

135.3 

-0.758 

0.272 

G2 

R3 

4452 

138.9 

-0.778 

0.272 

G3 

R1 

8935 

132.2 

-0.765 

0.264 

G3 

R2 

4780 

156.3 

-0.824 

0.264 

G3 

R3 

2470 

149.4 

-0.830 

0.264 

G4 

R1 

4571 

80.15 

-0.665 

0.303 

G4 

R2 

2297 

89.46 

-0.693 

0.303 

G4 

R3 

1174 

75.05 

-0.677 

0.303 
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Fig. 1. Wind tunnel. 
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b — O 




Fig. 2. Biplane grid configuration. 



Fig. 3. Model profile comparison. 
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R/U w dU/ds 



Figure 4. Elliptical leading edge models. 



Fig. 5. Stagnation point streamwise velocity gradient versus major to minor axis ratio. 



Fig. 6. Heat flux gage configuration. 
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Fig. 7. Components of total uncertainty in measured Frossling number. 
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s/R 


s/R 


a.) 1:1 ellipse 


c.) 2.25:1 ellipse 




s/R 


s/R 


b.) 1.5:1 ellipse 


d.) 3:1 ellipse 


Fig. 8. Low turbulence Frossling number distribution around leading edge compared to PARC-2D and Frossling solutions. 
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a.) 1:1 ellipse c.) 2.25:1 ellipse 




b.) 1.5:1 ellipse 


d.) 3:1 ellipse 


Fig. 9. Comparison of stagnation point heat transfer augmentation data and correlation. 
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Fig. 10. Distribution of averaged, normalized Frossling number downstream of stagnation point. 
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